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ABSTRACT  /Mox/muw  200 

The  primary  focus  of  the  investigation  was  to  conduct  a  chronological  evaluation  of  damage  initiation 
and  propagation  in  Electron  Beam  -  Physical  Vapor  Deposition  (EB-PVD)  thermal  barrier  coated 
(partially  stabilized  zirconia)  Rene  N5  single  crystal  superalloy.  Damage  initiation  and  failure  events 
due  to  thermal  cycling  were  investigated.  Early  crack  initiation  was  rationalized  through  interfacial 
asperity  induced  high  stresses  between  the  TBC  and  the  PtAl  bond  coat  .  Analytical  solution  was 
forwarded  for  transient  and  steady-state  residual  stresses  in  the  EB-PVD  system.  Oxide  scale 
development  as  a  function  of  time  was  predicted  using  Pick’s  Law.  The  deviation  of  experimental 
results  from  Pick’s  Law  was  attributed  to  microcrack  interaction  and  coalescence.  Oxide  scale  induced 
internal  pressure  on  crack  surfaces  was  shown  to  be  a  plausible  mechanism  for  microcrack  growth. 
Buckling  analysis  was  used  to  estimate  the  critical  size  of  delamination  necessary  (sixteen  times  the 
TBC  thickness)  for  spallation.  It  was  shown  that  for  various  thermal  cycles,  spallation  life  can  be 
estimated  based  on  critical  oxide  layer  thickness.  Thermal  wave  imaging  technique  was  used  to  track 
damage  condition  as  a  function  of  thermal  cycles.  It  was  found  that  saturation  of  thermal  wave 
amplitude  corresponds  to  spallation  life  of  the  sample.  This  work  provides  a  mechanism-based 
,  framework  for  life  prediction  in  TBC  systems. 
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INTRODUCTION 

Thermal  barrier  coatings  (TBCs)  provide  thermal  insulation  and  the  bond  coat  supports  a 
high  temperature  gradient  at  the  surface  of  high  temperature  alloy  substrates.  Application  of 
these  superalloy/TBC  systems  can  be  found  in  both  aerospace  and  land-based  gas  turbine 
engines.  In  automotive  applications,  the  piston  head  for  diesel  engine  is  coated  to  achieve  longer 
life  time  and  higher  performance  in  terms  of  fuel  reduction  and  power.  However,  these  coatings 
have  durability  problems,  due  to  the  material  and  thermal  mismatch  between  the  coating  and  the 
metallic  substrate.  Thermal  residual  stresses  develop  during  cool  down  from  processing 
temperatures  (Nusier,  and  Newaz,  1998;  Chang,  et  al.,  1987;  Williamson,  et  al.  1993;  Carapella 
et  al.,  1994;  Evans  et  al.,  1983).  Environmental  effects,  specifically  oxidation,  and 
thermomechanical  loading  create  additional  stresses  (Suo  1995;  Mikata,  and  Taya,  1985; 
Kardomateas,  1989;  Kardomateas,  1990;  Kardomateas,  and  Chung,  1994). 

These  stresses  may  initiate  microcracks  such  as  debonding  and  radial  cracks  and  can  have 
profound  effect  on  the  response  of  the  TBC  and  interfacial  damage  accumulation  and 
failure(Nusier,  and  Newaz,  1998;  Takeuchi  and  Kokini,  1994;  Homack  and  Kokini,  1988;  Evans 
and  Hutchinson,  1984;  Bottega  and  Maewal,1983;  Yin,  1985;  Chai,  1990;  Ryan  and  McCafferty, 
1995).  Their  understanding  is  essential  to  predict  the  behavior  of  the  coatings  and  their 
performance. 

Failure  processes  of  TBCs  have  received  considerable  attention  in  the  literature  as 
engineers  and  designers  have  made  efforts  to  improve  the  coating  manufacturing  process  and 
properties  to  improve  its  performance  in  a  variety  of  environments.  Partially  stabilized  zirconia 
(YSZ)  has  been  a  popular  TBC  for  numerous  applications.  The  bond  coat  is  often  the  oxidation 
barrier  to  the  substrate.  It  may  delay  but  can  not  prevent  oxidation  reaction.  Oxidation  of  the 
underlying  bond  coat  is  common  if  component  use  temperature  is  high.  Bond  coat  can  be  a 
number  of  material  systems  typically  denoted  as  MCrAlY  (M  stands  for  a  more  oxidation 
resistant  material  such  as  Nickel)  and  PtAl.  PtAl  diffusion  coating  has  been  developed  due  to 
demands  on  coatings  for  gas  turbine  engines  that  must  operate  at  significantly  higher 
temperatures  than  required  previously.  PtAl  is  also  cheaper  in  the  long  run  process.  The 
thermally  grown  oxide  (TGO)  layer  is  often  alumina  oxide,  AI2O3  which  grows  between  the 
TBC  and  the  bond  coat  layers.  Aluminum  is  consumed  from  the  bond  coat  to  form  AI2O3  as 
oxygen  diffusion  takes  place.  The  evolution  of  the  damage  process  including  oxidation  within 
the  TGO  layer  appears  to  control  the  overall  spallation  process.  However,  this  progressive 
damage  process  have  not  been  studied  adequately.  Rather,  researchers  have  proposed  a  number 
of  failure  scenarios  that  lead  to  spallation.  Because  of  inadequate  information,  either  a  physical 
or  a  phenomenological  basis  for  relating  spallation  life  to  growth  kinetics  within  the  TGO  layer 
or  damage  due  to  interlinking  microcracks  has  not  been  possible.  A  number  of  fundamental 
issues  remain  unanswered.  These  include  crack  initiation,  crack  propagation,  oxide  growth,  and 
total  damage  accumulation. 

There  are  three  most  plausible  scenarios  that  researchers  have  proposed  to  date  to  explain 
the  spallation  mechanism.  Evans  and  coworkers  (1983)  proposed  that  due  to  cool  down,  planar 
compressive  stresses  within  the  TBC  can  lead  to  buckling  provided  there  is  interfacial 
delamination  crack  already  present  between  the  TBC  and  the  substrate.  They  do  not  elaborate  on 
how  such  a  delamination  crack  may  occur  at  the  interface  in  the  first  place.  Our  recent 
calculations  with  the  current  EB-PVD  TBC  system  show  that  the  minimum  delamination  crack 
length  required  should  be  sixteen  times  the  TBC  thickness.  Buckling  of  the  TBC  layer  is 
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illustrated  in  Figure  1.  We  will  term  this  failure  model  as  Buckling  Model.  The  growth  of  the 
delamination  crack  is  possible  due  to  local  crack  tip  conditions  which  can  experience  both  shear 
and  out-of-plane  tension.  Another  model  for  spallation  is  due  to  another  author  Evans  (1989). 
The  basic  aspect  of  this  model  is  that  through-thickness  shear  crack  is  developed  in  the  TBC 
layer  under  compressive  conditions  in  the  TBC  as  shown  in  Figure  2.  This  model  can  be  termed 
as  Through-Thickness  Shear  Model.  Recently,  Suo  (1995)  proposed  the  Wrinkle  Model  which 
was  addressed  earlier  by  a  number  of  researchers  as  well  (Siegler,  1993;  Tolpygo  and  Grabke, 
1994).  According  to  this  failure  model,  there  is  void  formation  under  the  TBC  layer  and  folding 
effects  may  occur  that  lead  to  cracking  of  the  TBC  layer  at  various  locations. 

All  of  the  analysis  was  done  for  a  model  case  where  the  specimen  was  plasma-sprayed  in 
air  with  a  thin  zirconia-yttria  {ZrO^-iwt.VoY^O^)  layer  on  a  nickel-chromium-aluminum- 
zirconium  bond  coat.  The  properties  of  these  layers  are  given  in  Table  1. 

Characterize  progressive  damage  evolution  in  EB-PVD  TBC/superalloy  systems  leading 
to  spallation  and  to  develop  mechanics-based  models  to  rationalize  the  failure  events  are  the 
main  objectives  for  this  investigation. 

Firstly,  the  primary  focus  was  set  to  analytically  predict  the  transient  thermal  residual 
stress  distribution  at  the  various  interfaces  of  the  TBC/bond  coat/superalloy  system  and  compare 
the  results  with  finite  element  stress  analysis  results.  Evaluation  of  residual  stresses  is  treated 
until  it  reaches  a  steady-state  condition.  These  solutions  provide  an  insight  into  the  influence  of 
thermal  and  material  mismatch  stresses  at  the  interfaces  in  TBC  coated  superalloy  systems, 
particularly,  in  identifying  conditions  that  lead  to  interfacial  separation  due  to  thermal  loads. 

Secondly,  the  conditions  for  microcrack  initiation  due  to  thermal  residual  stresses  at 
wavy  interfaces  and  putting  an  empirical  equation  for  these  stresses  were  covered. 

Thirdly,  the  evaluation  of  the  conditions  necessary  for  growth  of  different  kinds  of 
interfacial  cracks  under  thermo-mechanical  loading  was  covered. 

Finally,  the  optical  microscopy,  scanning  electron  microscopy  and  thermal  wave  imaging 
techniques  were  used  to  characterize  the  interfacial  damage  of  a  two  layer  EB-PVD  TBC  on  a 
single  crystal  superalloy  (Rene’  N5)  substrate.  Pick’s  law  was  used  to  describe  the  TGO  buildup 
for  early  cycles. 

TRANSIENT  RESIDUAL  STRESSES  IN  THERMAL  BARRIER  COATINGS: 
ANALYTICAL  AND  NUMERICAL  RESULTS 

Elasticity  based  modeling  was  used  to  determine  the  transient  stresses  in  the  TBC,  the 
bond  coat ,  and  the  superalloy  substrate  with  specific  attention  to  the  interfaces.  For  the  steady 
state  case,  finite  element  modeling  was  imdertaken.  Closed  form  elasticity  solutions  correlated 
well  with  the  finite  element  results  for  the  steady  state  case.  The  highest  residual  stresses 
occurred  at  the  interface  between  the  bond  coat  and  the  TBC.  An  important  result  of  this 
investigation  was  that  the  TBC/bond  coat  interface  was  imder  biaxial  stress  field.  An  interesting 
result  was  that  the  residual  stresses  developed  in  the  substrate  are  higher  for  the  case  of  partly 
cooled  specimen  compared  to  the  fully  cooled  specimen  which  can  be  rationalized  due  to  the 
presence  of  higher  temperature  gradients  at  earlier  times  during  cool  down  from  processing 
temperature. 

The  orthogonal-expansion  technique  over  a  multilayer  region  was  used  in  order  to 
determine  the  temperature  distribution  in  the  concentric  cylinders.  The  first  three  roots  (coj^) 
were  found  to  be  dominant  in  prescribing  the  temperature  field.  All  other  roots  are  weak  and 
were  neglected.  The  axisymmetric  condition  was  applied  on  the  equilibrium  equations  along 
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with  stress-strain  relations  and  strain-displacement  relations,  so  the  equilibrium  equations  reduce 
to  two  ordinary  differential  equations.  The  transient  stress  field  was  obtained  hy  applying  the 
boundary  conditions  to  the  solution  obtained  from  those  ordinary  differential  equations. 
Transient  residual  stresses  have  been  evaluated  for  the  case  of  three  concentric  cylinder  by  finite 
element  approach  (steady  state  case  only)  as  well  as  using  elasticity  approach.  The  analytical 
model  is  described  in  detail  in  the  appendix  as  published  in  a  refereed  journal  recently  (Nusier, 
andNewaz,  1998). 

Figure  1  shows  the  variation  of  radial  stress  versus  radial  distance  at  different  time.  The 
peak  value  occurred  close  to  the  center  at  t=16.5  second.  Also  Figure  1  shows  that  the  radial 
stress  vanish  at  the  outer  surface.  Radial  stresses  are  not  critical  for  this  specimen  configuration 
because  of  their  low  values.  Note  that  radial  stress  range  is  0-25  MPa.  Figure  2  shows  the 
variation  of  tangential  stress  versus  radial  distance  at  different  time,  from  this  figure  one  can  see 
that  the  maximum  stress  occurred  in  the  bond  coat.  Figure  3  shows  the  variation  of  axial  stress 
versus  radial  distance  at  different  time,  same  as  in  Figure  2  the  maximum  axial  stress  occurred  in 
the  bond  coat.  These  figures  show  excellent  agreement  between  finite  element  solution  and 
elasticity  solution  for  the  steady  state  case.  The  bond  coat/TBC  interface  is  clearly  under  biaxial 
stress  field.  The  highest  tensile  stress  is  242  MPa  and  the  highest  compressive  stress  is  132  MPa. 
The  elasticity  solution  is  not  valid  close  to  the  free  edge  since  the  boundary  conditions  are  not 
satisfied  precisely  at  the  free  edge.  The  boundary  conditions  are  subscribed  in  average  form  at 
£iny  cross-section  which  may  not  be  precisely  valid  at  the  free  edge. 

CRACK  INITIATION  IN  THERMAL  BARRIER  COATINGS  DUE  TO  INTERFACE 
ASPERITY 

Microcracks  have  been  observed  to  initiate  near  the  bond  coat/TGO  interface  after  only  a 
few  thermal  cycles  in  electron  beam  -  physical  vapor  deposition  (EB-PVD)  TBC  systems.  The 
origin  of  these  microcracks  can  be  traced  to  the  buildup  of  thermal  stresses  which  may  magnify 
due  to  asperity  of  the  thermally  grown  oxide  (TGO)  layer  and  the  bond  coat.  The  TGO/bond  coat 
interface  was  modeled  as  rough  but  periodic  asperity.  Finite  element  modeling  was  undertaken  to 
determine  the  magnitude  of  stresses  at  a  rough  bond  coat/oxide  interface  (modeled  as  a  sine 
wave).  Large  normal  interface  stresses  arise  at  the  peaks  of  the  rough  interface,  while  large  shear 
stresses  arise  at  the  mean  line  of  the  rough  interface.  It  was  foxmd  that  these  residual  stresses  can 
initiate  microcracks  at  the  bond  coat/oxide  interface  that  can  result  in  debonding  and  can  also 
develop  cracks  within  the  TGO  layer.  High  wave  amplitude  and  small  wave  length  representing 
asperity  create  high  normal  and  shear  stresses  at  oxide  peak  and  at  the  mean  line  of  rough 
interface,  respectively.  In  particular,  the  condition  for  microcrack  initiation  was  investigated.  The 
results  show  that  asperity  can  result  in  stresses  high  enough  to  initiate  microcracks  at  the 
TGO/bond  coat  interface. 

The  interface  was  modeled  as  a  sine  wave  in  the  following  form 
Sir)  =  a,  sm(^)  (1) 

A- 

where  S(r)  is  the  interface  equation,  is  the  wave  amplitude,  r  is  the  radial  coordinate,  and  X  is 
the  wave  length.  The  stresses  of  interest  are  the  normal  and  shear  stresses.  Note  that  for  the  case 
of  flat  surface  at  the  interface,  the  normal,  and  shear  stresses  are  zero.  On  the  other  hand  the 
tangential  and  hoop  stresses  are  present,  and  they  are  very  large  in  the  oxide  layer  (4-5  GPa 
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compressive  stress).  Since  the  normal  and  shear  interface  stresses  are  acting  directly  on  the 
interface,  they  should  be  continuous  across  the  interface.  Tangential  and  hoop  stresses  act  within 
each  layer  on  either  side  of  the  interface,  and  will  therefore,  be  different  in  each  material. 

The  uncoated  superalloy  specimen  had  a  thickness  of  3.175  mm,  the  bond  coat  layer 
thickness  is  24.15  pm,  the  oxide  layer  thickness  is  24.15  pm,  the  ceramic  layer  (TBC)  thickness 
is  127  pm,  and  the  radius  of  the  specimen  is  12.7  mm.  The  specimen  was  cooled  down  from  a 

processing  temperature  of  1025  to  room  temperature  of  25  °C. 

A  parametric  study  was  carried  out  in  which  the  effect  on  interface  stresses  due  to 
variation  in  three  wave  amplitudes,  and  four  wave  lengths  representing  different  levels  of 
asperity  were  considered.  In  this  study  the  material  properties  were  considered  to  be  temperature 
dependent,  so  that  the  results  will  be  more  practical. 

For  the  case  of  2.413  pm  wave  amplitude  and  31.75  pm  wave  length,  the  peak  value  of 
interface  normal  stress  reaches  a  value  of  approximately  600  MPa  which  is  too  high  compared 
to  oxide  tensile  strength  (=215  MPa).  For  a  radial  length  of  a  bout  10  pm,  the  interface  normal 
stress  is  exceeding  the  tensile  strength  of  the  oxide  layer,  so  a  cracks  of  10  pm  will  be  created  at 
the  peaks  of  the  rough  interface.  Since  at  the  tips  of  these  cracks,  the  normal  stresses  are  still 
high,  one  can  expect  that  these  cracks  will  propagate  to  longer  ones.  On  the  other  hand,  the  shear 
stress  has  a  peak  value  at  the  mean  line  of  the  rough  interface,  for  a  case  similar  to  the  previous 
one,  cracks  of  15  pm  or  longer  will  be  created  at  the  mean  line  of  the  rough  interface,  since  the 
tensile  shear  strength  of  the  oxide/bond  interface  is  approximately  70  MPa.  Figures  4-5  show 
that  the  variation  of  normal  stresses  versus  wave  amplitude  and  the  inverse  of  the  wave  length  is 
close  to  be  linear. 

Based  on  the  previous  discussion  the  following  empirical  formulas  seems  to  be  a  good  fit 
for  normal  and  shear  stresses 
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where  s  is  the  thermal  strain,  v  is  Poisson’s  ratio,  E  is  the  modulus  of  elasticity,  and  subscripts 
o  and  b  corresponding  to  oxide  and  bond  coat  layers,  respectively.  Figure  6  shows  finite 
element  results  in  comparison  to  the  suggested  empirical  equations  for  normal  stresses.  Clearly 
the  empirical  equations  are  in  good  agreement  with  finite  element  results. 


ANALYSIS  OF  INTERFACIAL  CRACKS  IN  A  TBC/SUPERALLOY  SYSTEM  UNDER 
THERMAL  LOADING 

Effect  of  interfacial  microcracks  have  been  investigated  using  fracture  mechanics 
approach.  In,  particular,  J-integral  and  the  energy  release  rate  G,  for  both  Mode  I  and  Mode  II 
using  virtual  crack  extension  method  were  evaluated.  Two  types  of  specimens  were  studied.  The 
specimens  were  cooled  down  from  processing  temperature  of  1025  °C  to  25  ®C.  The  variation  of 
the  properties  as  function  of  temperature  were  used  for  analysis.  It  was  found  that  the  use  of 
temperature  dependent  properties  in  contrast  to  constant  properties  provide  significantly  different 
values  of  J-integral  and  G.  For  stepped-disc  specimen  with  edge  crack,  crack  growth  is  only  due 
to  Mode  II,  while  for  cylinder  specimen  with  internal  crack,  crack  growth  is  due  to  mixed-mode 
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loading.  An  important  implication  of  this  result  is  that  edge  delaminations  in  a  disk  specimen 
may  only  grow  due  to  Mode  II  conditions  imder  pure  thermal  loading.  Shear  fracture 
charaeteristics  of  interfacial  cracks  thus  become  important  in  the  failure  of  the  TBC. 

The  presence  of  circumferential  craek  between  the  ceramic  layer  and  the  bond  coat  layer 
was  analyzed  using  finite  element  method.  The  general  code  ABAQUS  was  used  for  these 
analysis.  The  energy  release  rate  G  was  estimated  by  using  the  virtual  crack  extension  method. 

Finite  element  method  was  used  in  order  to  determine  the  J-integral  value  and  the  energy 
release  rate.  Virtual  crack  extension  method  was  used  to  evaluate  the  energy  release  rate  for  both 
Mode  I  and  Mode  II.  For  cylinder  type  speeimen  the  imeoated  Nickel  based  superalloy  had  a 
radius  of  9.26  mm,  the  bond  coat  layer  thickness  is  0.184  mm,  the  TBC  thiekness  is  0.556  mm, 
and  the  length  of  the  specimen  is  100  mm.  For  stepped-disk  specimen,  the  disk  radius  is  40  mm, 
the  bond  coat  layer  thiekness  is  0.184  mm,  the  TBC  thickness  is  0.556  mm,  and  the  uncoated 
Nickel  based  superalloy  had  a  thickness  of  9.26  mm. 

Figure  7  shows  the  variation  of  J  and  G,  versus  crack  length  ratio  for  cylindrical  type 
specimen.  Clearly  both  of  them  have  a  very  elose  value,  also  J  and  G,  have  a  peak  value  higher 
than  the  steady  state  value  at  a  smaller  erack  length.  The  critical  value  of  the  energy  release  rate 
for  TBC  is  about  300-320  N.m/m^  based  on  fracture  toughness  of  pure  eeramic  material,  so 
xmder  applied  thermal  loading  the  crack  length  exceeds  four  times  the  TBC  thickness, 
instantaneous  craek  growth  leading  to  fracture  may  occur.  The  variation  of  J  and 
G,j{Gi  =G;;)  versus  crack  length  ratio  for  stepped-disk  specimen  is  shown  in  figure  8.  Both 
values  are  very  close,  and  for  a  very  short  erack  about  the  order  of  TBC  thickness,  the  J  and  Gj, 
has  a  significant  value  of  about  140  N.m/m^ .  Note  that  for  stepped-disk  speeimen  under  thermal 
loading.  Mode  I  is  essentially  absent  and  primarily  Mode  II  exists. 

ANALYSIS  OF  INTERFACIAL  CRACKS  IN  A  TBC/SUPERALLOY  SYSTEM  UNDER 
THERMO-MECHANICAL  LOADING 

Effeet  of  oxidation  growth  between  bond  coat  and  ceramic  layer  (TBC)  can  be  modeled 
as  volume  increase.  In  this  part  of  the  investigation,  we  represent  this  change  in  volume  as  an 
induced  pressure  across  the  crack  faces.  Mixed-mode  fracture  analysis  of  a  thin  circular 
delamination  in  an-axisymmetrically  multi  layer  circular  plate  is  developed.  Geometrical 
nonlinearity  is  included  in  the  analysis,  since  we  have  a  large  deflection  case.  The  elastic 
deformation  problem  of  a  circular  plate  subjected  to  a  clamped  boundary  condition  at  the  edge 
delamination,  an  out  of  plane  pressure  load,  and  a  compressive  stress  due  to  thermal  mismatch 
between  different  layers,  was  solved  numerically  using  a  Rayleigh-Ritz  method.  The  strain 
energy  release  rate  was  evaluated  by  means  of  the  path-independent  J-integral.  The  numerical 
results  of  this  problem  based  on  the  energy  method  were  verified  using  finite  element  method. 
Both  methods  correlate  well  in  predicting  the  energy  release  rate  for  Mode  I  and  Mode  II, 
deflection,  and  postbuckling  solutions.  The  energy  release  rate  G,  for  both  Mode  I  and  Mode  II 
cases  were  evaluated  using  the  virtual  crack  extension  method.  A  foxxr  layer  model  which  include 
the  TGO  layer  was  analyzed  using  finite  element  method,  two  cases  were  considered.  Case  one 
is  a  specimen  with  a  crack  at  middle  of  the  oxide  layer  (TGO)  while  case  2  is  a  specimen  with  an 
interface  crack  between  bond  coat  and  TGO  layers.  Also,  for  both  cases,  a  sensitivity  study  for 
the  effect  of  variation  of  materials  properties  (+25%)  were  undertaken  using  the  finite  element 
method. 

The  presence  of  cireumferential  craek  between  the  ceramic  layer  and  the  bond  coat  layer 
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was  analyzed  using  A  Rayleigh-Ritz  method  and  finite  element  method.  The  general  code 
ABAQUS  was  used  for  finite  element  analysis.  The  energy  release  rate  G  was  estimated  by 
using  the  virtual  crack  extension  method. 

Finite  element  method  was  used  in  order  to  determine  the  J-integral  value  and  the  energy 
release  rate.  Virtual  crack  extension  method  was  used  to  evaluate  the  energy  release  rate  for  both 
Mode  I  and  Mode  II  crack  growth.  For  stepped-disk  specimen,  the  disk  radius  is  12.7  mm,  the 
bond  coat  layer  thickness  is  0.04826  mm,  the  TBC  thickness  is  0.127  mm,  and  the  uncoated 
Nickel  based  superalloy  had  a  thickness  of  3.175  mm.  The  specimen  was  cooled  down  from  a 
processing  temperature  of  1025  ®C  to  a  temperature  of  25  “C . 

Figure  9  shows  the  variation  of  Gj  and  Gjj  versus  pressure  ratio  for  a  delamination 
radius  of  twenty  times  the  TBC  thickness.  This  figure  shows  good  agreement  between  finite 
element  method  results  based  on  VCEM  and  numerical  results.  Mode  I  is  dominant  in  this  case. 
Figure  10  shows  the  variation  of  J-integral  delamination  radius  from  numerical  analysis.  From 
this  figure  one  can  obtain  the  conditions  when  the  crack  will  propagate.  From  literature  the 

critical  energy  release  rate  is  varied  between  100-300  N.m/m^  for  this  system  configuration.  For 

a  value  of  300  J.m/m  ^ ,  and  a  delamination  radius  of  four  times  the  TBC  thickness,  the  crack  will 
propagate  at  induced  pressure  value  equivalent  to  500  atmosphere  ( «  50  MPa).  This  value  is  easy 
to  develop  due  to  volume  increase  because  of  oxide  layer  growth.  Also,  based  on  finite  element 
analysis  conducted  by  the  authors  for  this  system  with  a  wavy  interface,  the  results  shows  that  an 
axial  stress  of  1 85  MPa  can  be  developed  for  a  sine  wave  interface  with  amplitude  of  2.4  pm  and 
wave  length  of  127  pm. 


DAMAGE  ACCUMULATION  MECHANISMS  IN  THERMAL  BARRIER  COATINGS 
(EXPERIMENTAL  AND  ANALYTICAL  EVALUATION) 

Spallation  is  a  major  failure  condition  experienced  by  thermal  barrier  coatings  (TBCs) 
subjected  to  thermal  and  mechanical  loads.  Although  evidence  of  spallation  is  substantiated  and 
mechanistic  models  to  describe  the  failure  condition  is  prevalent  in  literature,  the  progressive 
nature  of  damage  evolution  leading  to  spallation  has  not  been  addressed  adequately.  In  this 
chapter,  discussion  on  the  investigation  of  damage  evolution  in  partially  stabilized  zirconia  TBC 
on  Nickel-based  single  crystal  superalloy,  Rene  N5.  Thermal  cycles  were  imposed  on  button 
and  rod  specimens  with  Electron  Beam  -  Physical  Vapor  Deposition  (EB-PVD)  TBC  coating. 
The  bond  coat  was  PtAl.  Different  thermal  cycles  profile  were  used  for  damage  evaluation 
purpose.  Progressive  damage  evolution  was  tracked  using  microscopy  eind  thermal  wave  imaging 
technique  on  samples  subjected  to  a  series  of  thermal  cycles.  Pick’s  law  was  used  to  describe  the 
thermally  grown  oxide  (TGO)  buildup  for  early  cycles.  It  was  found  that  Pick’s  law  can  be  used 
to  correlate  between  damage  thickness  and  oxide  thickness  for  different  thermal  cycles  profile. 
The  importance  of  this  correlation  is  that,  for  a  given  thermal  cycle  profile  and  number  of  cycle, 
the  oxide  thickness  can  be  found  anal5dically  using  Pick’s  law,  from  this  thickness  one  can 
determine  the  state  of  damage  (damage  thickness)  by  using  the  correlation  formula.  However,  at 
higher  number  of  thermal  cycles,  damage  in  the  form  of  microcracks  and  their  coalescence 
results  in  the  loss  of  integrity  of  the  TGO.  Thus,  both  oxidation  kinetics  and  damage  appears  to 
have  significant  roles  to  play  as  it  relates  to  spallation.  As  these  microcracks  coalesce  to  form 
major  delamination  cracks  or  interlayer  separation,  the  susceptibility  for  coating  buckling  is 
increased.  The  delamination  cracks  finally  consume  the  TGO  layer.  The  loss  of  TBC  integrity 
from  the  bond  coat  and  the  substrate  facilitates  its  buckling  during  cool  down  from  elevated 
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temperature.  Our  estimations  show  that  a  delamination  crack  length  of  about  twelve  to  sixteen 
times  the  TBC  thickness  is  needed  for  the  current  material  system  to  initiate  buckling. 
Progressive  microcrack  linking  is  a  possible  mechanism  to  develop  such  critical  delamination 
crack  lengths.  Physical  evidence  of  buckling  was  found  in  specimens  prior  to  complete 
spallation. 

EXPERIMENTAL  ASPECTS 

A  rapid  temperature  CM  furnace  (model  1610  BL  (c) )  was  used  for  thermal  cycling  tests. 
A  series  of  test  according  to  table  2  was  done  to  xmderstand  the  damage  accumulation 
mechanism  in  disk  top  type  specimens.  Each  cycle  has  a  heat  up  time  of  9  min.,  a  cooled  down 

time  of  10  min.,  and  a  of  200  "C.  After  each  test,  the  specimen  was  potted,  sectioned, 
grinded,  polished,  and  microscopy  pictures  were  taken,  respectively.  For  some  samples,  thermal 
wave  images  were  taken  in  to  quantitatively  identify  the  state  of  damage  (location  and  size). 
Microscopy  will  provide  some  information  on  the  damage  state.  Thermal  and  thermomechanical 
cycling  tests  on  cylindrical  type  specimens  were  conducted  by  using  the  CM  furnace  and  the 
hydraulic  MTS  machine  with  a  heating  control  unit  that  was  added,  respectively. 

THERMAL  CYCLING 

Thermal  cycling  tests  on  disk  top  and  cylindrical  type  samples  were  conducted  on  a  rapid 
temperature  CM  furnace  (model  1610  BL  (c)).  The  heating  zone  for  this  furnace  is  about  20  cm 
X  20  cm  X  25  cm.  To  insure  the  uniformity  of  the  temperature,  four  different  thermocouples 
(Type  K)  were  used,  two  of  them  located  at  the  bottom  of  the  furnace  and  they  are  very  close  to 
the  samples.  The  other  two  thermocouples  were  located  at  the  top  of  the  furnace.  A  digital 
recorder  (Honeywell)  were  used  to  record  the  temperature  variation.  It  should  be  mentioned  that, 
the  four  thermocouples  were  showing  a  temperature  difference  of  about  3-4  at  temperature  of 
1177  “C.  A  micristar  controller  was  used  to  set  the  parameters  used  to  define  a  thermal  cycle, 
such  as  heat  up  time,  hold  time,  cool  down  time,  minimum  temperature,  and  maximum 
temperature,  etc. 

THERMOMECHANICAL  CYCLING 

Thermomechanical  cycling  tests  conducted  on  cylindrical  type  specimen  were  under 
taken,  since  this  type  of  specimen  is  similar  to  turbine  blades.  The  specimen  will  be  imder 
compressive  load  generated  by  a  (55  Kip)  hydraulic  MTS  machine  fully  controlled  by  a  PC  (P5- 
60)  computer.  Heating  control  unit  was  added  to  support  the  required  temperature  variation  with 
time.  The  dimensions  of  cylindrical  type  specimen  used  are  5.9  mm  in  diameter  and  60  mm  in 
length. 

ANALYTICAL  ASPECTS  (Growth  of  TGO  Layer) 

The  development  of  the  diffusion  model  and  the  associated  governing  equations  is 
necessary  to  determine  the  time  dependent  oxide  layer  thickness.  Diffusion  is  approached  as  the 
penetration  of  one  substance,  such  as  oxygen,  into  a  porous  solid  matrix  of  a  different  substance. 


dc 

~dt 


=  (Dc,), 


(3) 
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Eqn.  3  is  Pick’s  second  law,  where  D  is  defined  as  the  diffusion  coefficient  {\im  Z^),  c  is  the 

oxygen  concentration  {atom /  |a.w^ ),  and  t  denote  time  dependence.  The  temperature  dependence 
is  incorporated  into  the  diffusion  coefficient  which  is  approximated  by  an  Arrhenius  type  relation 
ship 


D=D„^xp(^)  (4) 

where  Dq  is  the  diffusion  constant  ( /  s),Q  is  the  activation  energy  for  oxygen  diffusing  into 
metallic  substrate  (J/mol),  R  is  the  universal  gas  constant  (8.134  J/mol  K)  and  T  is  the 
temperature  (K).  (6.16) 

The  oxidation  growth  law  will  be  obtained  from  the  experiments  as  function  of  Df  and 
time,  where  is  the  effective  oxidation  constant  and  is  defined  as 

where  T(t)  is  the  temperature  which  can  vary  with  time. 

The  solution  of  the  partial  differential  Eqn.  3  for  a  one  dimensional  problem  is  obtained 
by  using  separation  of  variables  technique  in  the  following  form 

c(x,0  =  C,-(C,-C„)|;«,cosX,(i/i)c-‘’-""‘’  (6) 

n=0 


where, 

g  4sinX„ 
2X,„+sin2X,„ 


(7) 


The  eigenvalues  are  the  roots  for  the  following  equation, 

X„tan?i„-a,Z/Z)  =  0  (8) 

The  constants  used  for  analysis  in  this  investigation  were  as  follows  (Wright,  1996): 

Activation  energy  for  oxygen  diffusion:  Q  =  321  kJ/mol 
Diffusion  constant:  Do  =  3.85JflO'^  mVsec 
Universal  gas  constant:  R  =  8.314  J/mol  K 

The  diffusion  equation  for  along  a  cylinder,  and  for  the  case  were  the  variation  of 
concentration  with  the  angle  0  is  zero  (axisymmetric  case)  is  given  by 


dc  1  dc. 

—  =  D( — 7  + - ) 

dt  dr^  r  dr 


(9) 


The  solution  of  the  partial  differential  Eqn.  9  for  a  one  dimensional  problem  is  obtained 
by  using  separation  of  variables  technique  in  the  following  form 
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c(x,0  =  C„  -  2(C,  -  CjI;  /  -MX„r  / 

n=\  ^n^\y^n) 


(10) 


where  is  Bessel  function  of  the  first  kind  (order  n).  The  eigenvalues  are  the  roots  for  the 
following  equation, 


=  o  (11) 

Bond  coat  oxidation  resistance  has  been  clearly  linked  to  increased  durability  of  the  TBC. 
RESULTS  AND  DISCUSSION 

OPTICAL  MICROSCOPY  /  SCANNING  ELECTRON  MICROSCOPY  (SEM) 

The  thermal  cycle  profiles  used  for  conducting  all  tests  for  this  research  is  shown  in 
Figure  11.  Specimens  were  tested  for  various  thermal  cycles  and  examined  after  the  thermal 
exposures.  Specimens  were  evaluated  at  different  number  of  thermal  cycles  according  to  table  2. 
Complete  spallation  was  observed  in  these  specimens  between  175-185  cycles  (cycle  A),  435- 
455  cycles  (cycle  B),  and  100-205  cycles  (cycle  C).  The  thermally  cycled  samples  were  potted 
by  using  cold  mounting  epoxy  resin.  After  curing  for  one  day,  the  specimens  were  cut  by  using  a 
Buheler  Isomet  2000  diamond  cutter.  The  samples  were  then  ground  and  polished  by  using 
normal  metallurgical  procedvues.  The  polished  samples  were  observed  under  Olympus  BX  60 
optical  microscope.  A  few  samples  were  observed  imder  a  Hitachi  2000  Scanning  Electron 
Microscope  (SEM)  at  higher  magnifications  to  observe  the  nature  of  delamination  cracking 
within  the  TGO  layer. 

The  formation  of  an  oxide  layer  can  be  easily  traced  via  microscopy  at  the  interface 
between  the  TBC  layer  and  the  bond  coat.  Figure  12  shows  the  photomicrograph  of  a  thermal 
barrier  coated  sample  without  any  thermal  cycling.  As  expected,  no  interfacial  cracking  and 
thermally  grown  oxide  layer  were  evident  at  the  bond  coat/top  coat  (TBC)  interface.  Thermal 
cycling  of  samples  in  air  leads  to  the  formation  of  a  reaction  product  layer  and  interfacial 
cracking/damage  at  the  bond  coat/top  coat  interface.  This  oxidation  product  (TGO)  is  and 

has  been  confirmed  by  many  researchers.  The  TGO  appears  to  grow  thicker  with  higher  number 
of  thermal  cycles.  Also,  within  the  TGO  layer,  numerous  microcracks  and  void  like 
discontinuities  are  formed.  A  few  thermal  cycles  leads  to  the  formation  of  interfacial  cracking  at 
the  bond  coat/top  coat  interface  and  with  increasing  the  number  of  cycles  leads  to  the  formation 
of  both  interfacial  cracking  and  TGO  layer.  Figure  13  show  the  photomicrographs  of  the 
interfacial  cracking  at  the  bond  coat/top  coat  interface  taken  after  2  thermal  cycles  (cycle  type 
A).  Figures  14  through  17  show  the  chronological  evaluation  of  interfacial  damage  in  TBC  due 
to  thermal  cycling  at  25,  50,  110  and  175  (cycle  type  A).  The  thickness  of  the  TGO  layer  grows 
with  increasing  the  number  of  cycles. 

After  about  50  cycles  (cycle  type  A),  the  TGO  shows  separation  in  the  form  of 
delamination.  The  separation  between  the  TBC  and  the  bond  coat  takes  place  through  the  TGO 
layer  instead  of  delamination  at  any  preferred  interface.  The  separation  width  or  damage  width 
between  the  TBC  and  the  bond  coat  continues  to  increase  as  number  of  cycles  are  increased.  We 
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can  address  this  cycle  induced  change  by  accounting  for  damage  coalescence  within  the  TGO 
layer  and  increased  propensity  for  buckling  of  the  TBC  layer.  We  address  this  mechanism  next. 

Based  on  our  analysis  of  top  layer  buckling  on  a  circular  plate,  our  estimates  show  that 
buckling  is  possible  only  if  the  delamination  crack  size  is  twelve  to  sixteen  times  or  higher  than 
the  TBC  thickness.  The  evidence  of  buckling  in  a  specimen  cycled  to  175  (cycle  type  A)  cycles 
where  the  separation  vddth  for  delamination  is  considerably  higher  at  the  middle  of  the  specimen 
than  near  the  ends  points  out  that  the  delamination  crack  through  the  TGO  layer  may  reach  such 
a  critical  length  necessary  to  induce  buckling  .  The  evolution  of  damage  in  the  form  of 
microcracks  and  voids  within  this  layer  can  interact  to  form  such  delamination  cracks. 

For  only  oxidation  driven  degradation,  TGO  layer  growth  can  be  described  by  Pick’s  law. 
The  predicted  TGO  layer  growth  is  shown  in  Figures  18  and  19  for  cycle  A  &  B  respectively. 
The  TGO  layer  growth  in  addition  to  separation  thickness  is  termed  in  Figures  18  and  19  as 
“damage  thickness”.  Damage  thickness  in  this  study  is  referred  to  any  change  in  specimen 
structure,  these  changes  can  be  oxide,  separation,  crack,  voids,  ...etc.  Based  on  that  when  we 
have  breakdown  of  the  TGO  layer,  we  refer  to  the  difference  between  damage  thickness  and 
oxide  thickness  as  separation.  The  experimental  data  which  is  an  average  of  at  least  four  optical 
micrographs  for  each  point  along  the  length  of  the  TGO  layer  is  shown  in  the  same  figure.  A 
computer  code  (written  in  C** )  was  developed  to  evaluate  the  damage  area  from  which  the 
damage  thickness  can  be  evaluated.  The  input  for  this  program  is  the  output  of  a  scanned  optical 
micrograph  in  TIF  format. 

The  deviation  of  experimental  results  from  the  predicted  TGO  layer  thickness  occurs  at 
about  50  thermal  cycles  (cycle  A)  and  at  130  thermal  cycles  (cycle  B).  According  to  Pick’s  law, 
these  two  cycles  have  approximately  the  same  oxide  thickness.  The  significance  of  this  result  is 
that  below  50  thermal  cycles  of  the  type  chosen  in  this  investigation,  the  TGO  layer  growth  is 
oxidation  driven. 

According  to  Eqn.  5,  for  any  thermal  cycle  profile,  an  effective  diffusion  constant  can  be 
evaluated  by  numerical  integration  (Simpson’s  rule  was  used).  After  the  effective  diffusion 
constant  has  been  determined,  the  oxidation  thickness  according  to  Eqn.  6  (button  samples)  and 
Eqn.  10  (cylindrical  specimen)  can  be  evaluated.  A  very  interesting  result  were  found,  that  is  the 
damage  thickness  (experimental  measurement)  is  equal  to  oxide  thickness  (Pick’s  law)  up  to  a 
critical  oxide  thickness  (=6.2  pm),  beyond  that  thickness,  the  daimage  thickness  start  to  get 
bigger  than  the  oxide  thickness  (separation  thickness  is  the  difference  between  damage  thickness 
and  oxide  thickness),  in  another  word  separation  start  to  occurred.  At  another  critical  oxide 
thickness  (si  3  pm)  the  sample  spall  off  This  observation  has  been  proven  for  different  thermal 
cycles  profiles,  different  maximum  temperature  and  different  holding  time  were  used  (cycles 
type  A,  B,  C  and  D).  Figure  20  shows  the  variation  in  total  damage  thickness  versus  oxide 
thickness  for  four  different  thermal  cycles.  This  figure  clearly  shows  that  four  different  thermal 
cycle  profile  follow  same  trend  to  failure.  An  important  outcome  from  this  figure,  is  that,  for  any 
thermal  cycle  profile  and  at  any  time  during  thermal  exposure,  the  oxide  thickness  can  be 
calculated  very  easy  and  straight  forward  by  the  use  of  Pick’s  law,  knowing  the  oxide  thickness 
will  lead  to  identify  the  damage  state  (Figure  20).  To  develop  Figure  20,  the  maximum 
temperature  were  varied  between  1177-1132  °C  (this  is  the  recommended  application 
temperature  range  for  this  kind  of  samples  by  GE),  and  the  holding  time  were  varied  between  45- 
600  minutes.  As  will  be  discussed  later  on,  the  holding  time  has  no  effect  on  damage  size  as  long 
as  the  sum  of  holding  time  is  fixed.  The  maximum  temperature  has  a  major  effect  on  life  time, 
for  example,  reducing  the  maximum  temperature  from  1177  ®C  to  1135  ®C  will  increase  the 
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sample  life  time  by  a  factor  of  2.4.  This  observation  is  illustrated  in  Figure  21,  this  figure  has 
been  developed  based  on  Pick’s  law  and  the  critical  oxide  thickness  to  failure  which  is  13  |a.m.  It 
is  very  clear  that  time  to  failure  is  highly  dependent  on  the  maximum  temperature.  A  very 
important  result  that  is  obtained  from  this  research  is  that,  the  oxide  thickness  evaluated  by 
Pick’s  law  can  be  still  be  used  as  a  correlation  parameter  to  define  the  state  of  damage  for  any 
thermal  cycle  profile  and  at  any  time  during  thermal  exposure. 

For  the  case  of  cylindrical  type  specimens,  two  specimens  were  run  under  80  MPa 
compressive  stress  and  thermal  cycle  type  A.  The  first  one  failed  at  103  cycle  and  the  second 
sample  failed  at  112  cycle.  Under  the  same  conditions  with  zero  compressive  stress  a  third 
sample  failed  after  108  thermal  cycle.  The  compressive  stress  shows  no  effect  on  life  due  to  the 
fact  that  only  a  small  part  of  the  compressive  stress  was  carried  by  the  TBC  (=13  MPa).  Based 
on  these  three  samples  and  Pick’s  law  (Eqn.  10),  a  critical  oxide  thickness  to  failure  was 
determined  as  10.5  pm. 

Figure  22a  shows  the  thermal  wave  image  of  four  samples  at  0,  25,  50  and  100  thermal 
cycling  (cycle  A)  respectively  from  left  to  right.  The  corresponding  surface  temperature  profile 
plots  shown  on  the  right  hand  side  of  the  figure  represents  the  relative  amplitude  plotted  as 
function  of  position  or  distance  across  the  center  of  the  samples.  The  plot  indicates  that  the 
thermal  wave  signal  amplitude  increases  100%  from  0  thermal  cycling  to  25  thermal  cycling, 
while  the  change  in  the  thermal  wave  signal  amplitude  from  25  thermal  cycling  to  50  thermal 
cycling  is  only  10%  for  each  case.  Figure  22b  shows  the  thermal  wave  image  of  a  second  set  of 
samples  at  0,  150,  175  thermal  cycling,  respectively,  from  left  to  right.  The  corresponding 
surface  temperature  profile  shows  that  the  thermal  wave  signal  amplitude  increases  150%  from  0 
thermal  cycling  to  150  thermal  cycling,  however,  the  thermal  wave  signal  amplitudes  of  the  150 
thermal  cycling  and  175  thermal  cycling  samples  are  identical.  With  greater  number  of  thermal 
cycles,  the  damage  thickness  at  the  bond  coat/top  coat  interface  are  also  increased.  The  damage 
thickness  may  also  change  depending  on  the  number  of  thermal  cycles.  This  damage  thickness 
may  be  responsible  for  the  thermal  wave  signal  amplitude  to  vary. 

SUMMARY  AND  CONCLUSIONS 

•  Analysis  of  concentric  cylinder  model  representing  TBC/bond  coat/superalloy  system  showed 
that  elasticity  solution  and  numerical  results  from  finite  element  stress  analyses  (steady  state  case 
only)  correlate  well  in  predicting  thermal  residual  stresses  in  the  various  layers  and  at  the 
interfaces.  The  highest  stresses  were  found  to  occur  at  the  bond  coat/TBC  interface.  At  this 
interface,  results  indicate  a  biaxial  stress  field  due  to  high  axial  and  tangential  stress;  the  radial 
stress  component  was  relatively  small.  The  residual  stresses  developed  in  the  substrate  are  higher 
for  the  case  of  partly  cooled  specimen  compared  to  the  fully  cooled  specimen  because  the 
temperature  field  has  a  high  gradient  at  earlier  stages  compared  to  later  stages. 

•  Finite  element  modeling  was  imdertaken  to  determine  the  variation  in  residual  stresses  at  a 
rough  bond  coat/oxide  interface  (modeled  as  a  sine  wave).  Large  normal  interface  stresses  arise 
at  the  peaks  of  the  rough  interface,  while  large  shear  stresses  arise  at  the  mean  line  of  the  rough 
interface.  High  wave  amplitude  and  small  wave  length  signifying  different  conditions  of  asperity 
create  a  high  normal  and  shear  stresses  at  oxide  peak  and  at  the  mean  line  of  rough  interface, 
respeetively.  In  particular,  the  condition  for  microcrack  initiation  was  determined.  The  results 
show  that  asperity  can  result  in  stresses  high  enough  to  initiate  microcracks.  Suggested  empirical 
equations  for  normal  and  shear  stresses  are  in  good  agreement  with  finite  element  results. 
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•  Effect  of  interfacial  microcracks  have  been  investigated  using  fracture  mechanics  approach 
under  pure  thermal  loading.  The  use  of  temperature  dependent  properties  in  contrast  to  constant 
properties  provide  significantly  different  values  of  J-integral  and  G  values.  Edge  delaminations 
in  a  stepped-disk  specimen  may  only  grow  due  to  Mode  II  conditions  under  pure  thermal 
loading.  Cylindrical  type  specimen  with  a  circumferential  crack  has  mixed  mode  conditions 
under  pure  thermal  loading.  For  a  stepped-disk  specimen  with  a  very  short  edge  crack,  the 
energy  release  rate  and  J-integral  values  are  significantly  high.  This  is  expected,  because  of  the 
stress  singularity  at  the  bimaterial  interface  at  the  free  surface. 

•  Effect  of  interfacial  microcracks  was  investigated  using  fracture  mechanics  approach  under 
thermo-mechanical  loading.  Thermal  load  was  used  in  conjimction  with  internal  surface  induced 
pressure  for  interfacial  cracks.  Stepped-disk  type  specimen  with  a  central  crack  has  mixed  mode 
conditions  under  thermo-mechanical  loading.  The  maximum  deflection  and  total  energy  release 
rate  evaluated  by  the  finite  element  method  and  the  energy  method  are  in  very  good  agreement. 

For  a  TBC  toughness  value  of  300  J.m/m  ^ ,  and  a  delamination  radius  of  four  times  the  TBC 
thickness,  the  crack  will  propagate  at  induced  pressure  value  equivalent  to  500  atmosphere  (»  50 
MPa).  The  results  clearly  indicate  that  small  internal  pressure  due  to  oxidation  induced  volume 
change  may  create  the  necessary  conditions  for  crack  growth  dining  thermal  cycling.  This  value 
may  easily  develop  due  to  volume  increase  because  of  oxide  layer  growth  or  due  to  the  presence 
of  a  wavy  interface. 

•  Progressive  damage  evolution  was  tracked  using  microscopy  and  thermal  wave  imaging 
technique  on  samples  subjected  to  a  series  of  thermal  cycles.  It  was  found  that  Pick’s  law  can  be 
used  to  correlate  between  damage  thickness  and  oxide  thickness  for  different  thermal  cycles 
profile.  The  maximum  temperature  has  a  major  effect  on  lifetime,  for  example,  reducing  the 
maximum  temperature  from  1177  C°  to  1 135  C°  will  increase  the  sample  lifetime  by  a  factor  of 
2.4.  If  the  sum  of  the  holding  times  is  kept  constant,  then  it  leads  to  similar  damage  state,  that  is, 
isothermal  and  thermal  cycling  have  similar  damage  states  provided  that  the  peak  temperature  is 
the  same. 

•  Based  on  the  investigation  imdertaken  in  the  study,  initiation  of  microcracks  at  the  TGO/Bond 
coat  interface  can  be  related  to  asperity  induced  high  stresses.  Coalescence  of  these  microcracks 
can  be  related  to  their  growth  influenced  by  oxide  growth  which  may  produce  crack-surface 
pressure  under  thermal  cycling.  Finally,  the  coalesced  microcracks  can  form  a  major 
delamination  crack  at  the  TGO/Bond  coat  interface.  For  large  delamination  cracks  (over  12  times 
TBC  thickness),  in-plane  compressive  stress  in  the  TBC  layer  due  to  the  cool-down  part  of  the 
thermal  cycle  can  lead  to  spallation  of  the  TBC  layer.  Thus,  the  progressive  nature  of  damage 
initiation  and  accumulation  leading  to  TBC  failure  are  forwarded  in  this  study  and  are 
rationalized  through  mechanics  analyses. 
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TABLE  1 


0 

Material  properties  at  22,  566  and  1149  C 
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Table  2 

Test  matrix  for  disk  top  type  specimens 


Ns 

Nc 

KJc 

H.T 

min 

C.T 

min 

Hd.T 

min 

Spall 

off 

Sect. 

T.W.I 

1 

2 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

2 

18 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

3 

25 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

4 

42 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

5 

50 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

6 

75 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

7 

100 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

8 

110 

200 

1177 

9 

10 

45 

No 

Yes 

No 

9 

150 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

10 

175 

200 

1177 

9 

10 

45 

No 

Yes 

Yes 

11 

175 

200 

1177 

9 

10 

45 

Yes 

No 

No 

12 

185 

200 

1177 

9 

10 

45 

Yes 

No 

No 

13 

25 

200 

1132 

9 

10 

45 

No 

Yes 

No 

14 

50 

200 

1132 

9 

10 

45 

No 

Yes 

No 

15 

100 

200 

1132 

9 

10 

45 

No 

Yes 

No 

16 

150 

200 

1132 

9 

10 

45 

No 

Yes 

No 

17 

200 

200 

1132 

10 

45 

No 

Yes 

No 

18 

250 

200 

1132 

9 

10 

45 

No 

Yes 

No 

19 

300 

200 

9 

10 

45 

No 

Yes 

No 

20 

390 

200 

1132 

9 

10 

45 

No 

Yes 

No 

21 

435 

200 

1132 

9 

10 

45 

Yes 

No 

No 

22 

455 

200 

1132 

9 

10 

45 

Yes 

No 

No 

23 

40 

200 

1172 

9 

10 

60 

No 

Yes 

No 

24 

100 

200 

1172 

9 

10 

45 

Yes 

No 

No 

139 

200 

1172 

9 

10 

45 

No 

Yes 

No 

205 

200 

1172 

9 

10 

45 

Yes 

No 

No 

27 

50 

200 

9 

10 

60 

No 

Yes 

Yes 

28 

50 

200 

1172 

9 

10 

60 

No 

Yes 

Yes 

29 

5 

200 

1172 

9 

10 

600 

No 

Yes 

Yes 

30 

5 

200 

1172 

9 

10 

600 

No 

Yes 

Yes 

Where,  is  sample  number,  is  cycle  number,  H.T  is  heat  up  time,  C.T  is  cool  down  time, 
Hd.T  is  hold  time.  Sect,  stand  for  wither  the  specimen  has  been  sectioned  or  not,  and  T.W.I  is  the 
abbreviation  for  thermal  wave  imaging  technique. 
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Figure  1.  Radial  stress  versus  radial  distance  at  different  time 
for  the  TBC/bond  coat/superalloy  system 
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Figure  4.  Maxmium  interface  normal  stress  versus  wave 
amplitude. 


Maxmium  Interface  Normal  Stress,  MPa 


-  aQ=2.413  ^im 

- aQ=1.2065|am 

- ao=0.60325  |xm 


T  -1 

X  ,  i^m 


Figure  5.  Maxmium  interface  normal  stress  versus  the 
inverse  of  the  wave  length. 


Figure  6.  Interface  normal  stress  versus  radial  position. 
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Figure  7.  Variation  of  J  and  G,  versus  crack  length 
ratio  for  cylindrical  specimen 
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Cycle 

"^max 
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A 

1177 
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9 

54 

64 

B 

1132 

200 

9 

54 

64 

C 

1172 

200 

9 

69 

79 

D 

1172 

200 

9 

609 

619 

Figure  1 1 .  Thermal  cycling  profiles 
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Figure  12  shows  an  optical  micrograph  of  an  untested  thermal  barrier  coated  sample.  Thermally 
grown  oxide  (TGO)  layer  and/or  interfacial  seperation  were  absent  at  the  bond  coat/top  coat 
(TBC)  interface. 


Figure  13  shows  an  optical  micrograph  of  a  tested  thermal  barrier  coated  sample.  Thermally 
grown  oxide  (TGO)  layer  and/or  interfacial  seperation  were  seen  at  the  bond  coat/top  coat  (TBC) 
interface.  The  sample  was  thermally  cycled  to  2  cycle  (type  A). 
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Figure  14  shows  an  optical  micrograph  of  a  tested  thermal  barrier  coated  sample.  Thermally 
grown  oxide  (TGO)  layer  and/or  interfacial  seperation  were  seen  at  the  bond  coat/top  coat  (TBC) 
interface.  The  sample  was  thermally  cycled  to  25  cycle  (type  A). 
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Figure  16  shows  an  optical  micrograph  of  a  tested  thermal  barrier  coated  sample.  Thermally 
grown  oxide  (TGO)  layer  and/or  interfacial  seperation  were  seen  at  the  bond  coat/top  coat  (TBC) 
interface.  The  sample  was  thermally  cycled  to  1 10  cycle  (type  A). 


Figure  17  shows  a  middle  optical  micrograph  of  a  tested  thermal  barrier  coated  sample. 
Thermally  grown  oxide  (TGO)  layer  and/or  interfacial  seperation  were  seen  at  the  bond  coat/top 
coat  (TBC)  interface.  The  sample  was  thermally  cycled  to  175  cycle  (type  A). 


thickness,  tim 
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Figure  21 .  Variation  of  versus  time  to  failure  based  on 
Pick's  law  and  critical  oxide  thickness  (13  um). 
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Comparison  of  thermal  wave  images  at  various  thermal  cycling  and 
corresponding  surface  temperature  profile  across  the  specimens 


OCys  150  Cys  175  Cys 


Figure  22.  The  thermal  wave  image  of  six  samples  at  0, 25,  50,  100,  150  and  175  thermal 
cycles  (cycle  type  A)  and  their  corresponding  surface  temperature  profiles  across  the 
specimens. 
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Thermal  barrier  coatings  (TBCs)  provide  thermal  insulation  to  high-temperature 
superalloys.  Residual  stresses  develop  in  TBCs  during  cool-down  from  processing 
temperatures  due  to  the  thermal  expansion  mismatch  between  the  different  layers 
(substrate,  bond  coat,  and  the  ceramic  TBC).  These  residual  stresses  can  initiate 
microcracks  at  the  bond  coat/TBC  interface  which  can  lead  to  debonding  at  the 
bond  coat/TBC  interface.  Elasticity -based  modeling  was  used  to  determine  the  tran¬ 
sient  stresses  in  the  TBC,  bond  coat,  and  the  superalloy  substrate  with  specific 
attention  to  the  interfaces.  For  the  steady-state  case,  finite  element  modeling  was 
undertaken  as  well.  Closed-form  elasticity  solutions  correlated  well  with  the  finite 
element  results  for  the  steady-state  case.  The  highest  residual  stresses  occurred  at 
the  interface  between  the  bond  coat  and  the  TBC.  An  important  result  of  this  investiga¬ 
tion  was  that  the  TBC/bond  coat  interface  was  under  biaxial  stress  field.  An  important 
result  was  that  the  residual  stresses  developed  in  the  substrate  are  higher  for  the 
case  of  partly  cooled  specimen  compared  to  the  fully  cooled  specimen  which  can  be 
rationalized  due  to  the  presence  of  higher  temperature  gradients  at  earlier  times 
during  cool-down  from  processing  temperature. 


Introduction 

Thermal  barrier  coatings  (TBCs)  provide  thermal  insulation 
and  the  bond  coat  provides  oxidation  resistance  at  high  tempera¬ 
tures.  Ceramic  matrix  composites  (CMCs)  can  be  used  for  high- 
temperature  application,  whereas  metallic  substrate  such  as  su¬ 
peralloys  cannot  perform  at  this  high  temperature  without  this 
ceramic  layer,  i.e.,  the  TBC.  Plasma-sprayed  zirconia-yttria  ce¬ 
ramic  layer  with  a  nickel-chromium-aluminum-yttrium  bond 
coat  on  a  substrate  made  of  nickel-based  superalloy  (Chang  et 
al,  1987)  is  a  common  superalloy /TBC  system.  Application 
of  these  superalloys  can  be  found  in  aircraft  engine  (advanced 
gas  turbines)  parts  such  as  ducts,  vanes,  and  blades,  etc.  In 
automotive  applications,  the  piston  head  for  diesel  engine  is 
coated  to  achieve  longer  life  time  and  higher  performance  in 
terms  of  fuel  reduction  and  power.  However,  these  coatings 
have  durability  problems,  due  to  the  material  and  thermal  mis¬ 
match  between  the  coating  and  the  metallic  substrate.  Thermal 
residual  stresses  develop  during  cool  down  from  processing 
temperatures.  These  residual  stresses  may  initiate  microcracks 
such  as  debonding  and  radial  cracks  and  can  have  a  profound 
effect  on  the  response  of  the  TBC  and  interfacial  damage  accu¬ 
mulation  and  failure.  Their  understanding  is  essential  to  predict 
the  behavior  of  the  coatings  and  their  performance. 

A  quantitative  determination  of  the  stress  states  in  a  model 
thermal  barrier  coating  as  it  cools  down  from  the  processing 
temperature  (stress-free)  to  another  lower  temperature  was 
studied  by  Chang,  et  al.,  (1987).  Considering  a  three-layer 
concentric  cylinder  model  and  a  plane-strain  condition,  they 
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clearly  illustrated  the  importance  of  the  coefficient  of  thermal 
expansion  of  the  bond  coat  and  the  Young’s  modulus  of  the 
TBC.  The  finite  element  approach  was  used  by  these  authors. 
The  experimental  method  and  finite  element  approach  was  used 
by  Holmes  (1988)  to  study  thermal  fatigue  of  coated  super¬ 
alloys  of  a  stepped-disk  specimen.  The  results  obtained  for  the 
various  stress-strain  histories  clearly  illustrate  that  both  the  com¬ 
pressive  strains  encountered  on  specimen  heatup  and  the  tensile 
strains  encountered  on  cool-down  can  have  a  profound  effect 
on  useful  coating  life.  An  elastic-plastic  finite  element  method 
was  formulated  to  study  residual  stresses  developed  at  graded 
ceramic-metal  interfaces  during  cool-down  by  Williamson  et 
al.,  ( 1993).  The  results  demonstrate  the  importance  of  account¬ 
ing  for  plasticity  when  comparing  graded  and  nongraded  inter¬ 
faces.  TTiey  studied  the  effects  of  different  interlayer  thicknesses 
and  nonlinear  composition  profiles  on  strain  and  stress  distribu¬ 
tions  established  during  cool  down  (Drake  et  al.,  1993).  Mikata 
and  Taya  ( 1985 )  studied  the  stress  field  in  a  coated  continuous 
fiber  composite  subjected  to  thermomechanical  loading,  a  four- 
layer  concentric  circular  model  was  used.  The  displacement 
approach  was  used  by  Kardomateas  (1989)  to  study  the  tran¬ 
sient  thermal  stress  problem  for  a  hollow  circular  cylinder.  He 
presented  the  variation  of  stresses  as  function  of  time  and  thick¬ 
ness.  The  effects  of  circumferential  variation  of  the  heat  transfer 
coefficient  and  gas  radiation  in  three  concentric  cylinders  were 
studied  by  Chung  et  al.  (1985).  They  determined  the  transient 
two-dimensional  temperature  distribution  by  using  the  finite 
element  method.  Transient  thermal  stresses  in  fully  and  partly 
cooled  circular  rings  were  investigated  by  Iwaki  (1986).  He 
showed  that  the  thermal  stress  in  the  fully  cooled  circular  ring 
is  not  always  greater  than  that  in  the  partly  cooled  circular  ring. 
Thermal  fracture  of  multilayer  ceramic  thermal  barrier  coatings 
was  studied  by  Takenchi  and  Kokini  (1994).  Controlled  experi¬ 
ments  along  with  finite  element  modeling  were  used.  Their  work 
showed  that  surface  cracks  are  formed  as  a  result  of  tensile 
stresses  that  are  created  following  stress  relaxation  in  the  TBC 
at  steady-state  high  temperatures.  The  effect  of  a  transient  ther- 
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mal  load  on  a  coating  which  is  bonded  to  a  cylindrical  substrate 
was  studied  by  Kokini  and  Homack  (1988).  The  finite  element 
method  was  used  to  obtain  a  solution  for  a  circumferential  edge 
crack  normal  to  the  coating.  Their  analysis  showed  that  smaller 
heat  transfer  rates  result  in  smaller  stress  intensity  factors.  The 
initial  phase  of  transient  thermal  stresses  due  to  general  bound¬ 
ary  thermal  loads  in  an  orthotropic  hollow  cylinder  was  studied 
by  Kardomateas  ( 1990).  Kardomateas  and  Chung  ( 1994)  stud¬ 
ied  the  case  of  orthotropic  cylindrical  shell  loaded  by  an  external 
pressure  under  imposed  constant  moisture  concentrations  on 
the  inner  and  outer  surfaces.  They  obtained  an  exact  elasticity 
solution  for  stresses  and  displacements.  The  thermomechanical 
response  of  ceramically  coated  metal  parts  in  elevated  thermal 
environments  were  investigated  by  Padovan  et  al.,  (1984, 
1985).  A  finite  element  algorithm  was  used  in  their  analysis. 
The  thermomechanical  oxidation  response  of  thermal  barrier 
coatings  was  studied  by  Padovan  et  al.,  (1986).  They  solved 
the  field  equations  with  the  associated  boundary  conditions  nu¬ 
merically. 

A  thermal  barrier  coated  high-temperature  substrate  in  the 
form  of  a  cylinder  is  a  multilayer  system  with  a  number  of 
interfaces.  The  substrate  is  first  coated  with  a  bond  coat  whose 
primary  function  is  to  provide  oxidation  protection  and  to  en¬ 
hance  bonding  of  the  thermal  barrier  coating  (TBC)  to  the 
overall  system.  The  TBC  layer,  the  bond  coat,  and  the  substrate 
have  different  thermal  expansion  coefficients.  As  a  result  due 
to  thermal  loading  and  related  excursion,  each  layer  may  experi¬ 
ence  residual  thermal  stresses.  Prior  to  reaching  a  steady-state 
situation,  the  thermal  residual  stresses  can  vary  significantly 
over  time,  particularly,  at  early  times,  say  within  the  first  100 
seconds  of  cooling.  Therefore,  the  problem  of  transient  residual 
stresses  can  become  quite  important  and  can  affect  the  bond 
coat/TBC  and  the  substrate /bond  coat  interfaces.  Numerical 
solution  techniques  are  often  utilized  for  such  problems  pro¬ 
vided  that  the  geometry  of  the  multilayer  system  is  amenable 
to  idealization.  However,  for  TBC  systems  there  no  closed-form 
solutions  available  to  determine  the  overall  temperature  field 
and  subsequently  to  obtain  the  transient  residual  stress  field. 
In  the  current  work,  a  concentric  cylinder  model  is  taken  for 
representating  a  realistic  superalloy  substrate  coated  with  bond 
coat  and  TBC  and  to  obtain  an  elasticity  solution  that  describes 
the  transient  residual  stress  fields. 

In  our  investigation,  the  primary  focus  was  to  analytically 
predict  the  transient  thermal  residual  stress  distribution  at  the 
various  interfaces  of  the  TBC/bond  coat/ superalloy  system  and 
compare  the  results  with  finite  element  stress  analysis  results. 
Evaluation  of  residual  stresses  is  treated  until  it  reaches  a  steady- 
state  condition.  These  solutions  provide  an  insight  into  the  in¬ 
fluence  of  thermal  and  material  mismatch  stresses  at  the  inter¬ 
faces  in  TBC  coated  superalloy  systems,  particularly  in  identi¬ 
fying  conditions  that  lead  to  interfacial  separation  due  to  thermal 
loads. 


Modeling  Aspects 

For  the  case  of  N  concentric  cylinders,  it  is  possible  to  get 
an  exact  elasticity  solution  far  away  from  the  edges.  This  pro¬ 
vides  a  better  understanding  of  the  type  of  stresses  that  exist, 
and  also  we  will  be  able  to  verify  these  stresses  with  the  finite 
element  solution.  The  elasticity  solution  is  developed  as  follows. 

Transient  Temperature  Solution.  Consider  a  one-dimen¬ 
sional  multiphase  region  consisting  of  N  concentric  cylinders, 
where  the  outermost  one  is  layer  number  one  and  the  innermost 
one  is  solid  and  has  a  layer  number  N .  Continuity  of  temperature 
and  heat  flux  was  assumed  to  be  across  the  interfaces  (perfect 
thermal  contact  between  the  layers  at  the  interfaces ).  The  gov¬ 
erning  equation  for  the  transient  temperature  field  is 


d%  1  OTi  ^  J^dJ^ 
dr^  r  dr  a,  dt 

/  =  1,  2,  .  .  .  ,  r,+,  <  r  <  r,  f  >  0.  ( 1 ) 

The  boundary  conditions  are 

r,.,(r,>,,r)  =  r,(r,^,,r)  /-1, 2 . iV-1  t>0  (2) 

,  dTi^iin^x,  t)  ,  dTiin^ut) 

- Vr - 

1,2,  I  t>0  (3) 

hT,(rut)  +  =  0  t  >  0.  (4) 

dr 

The  initial  condition  is 


7^(r,  0)-ro  r,+  i  <  r  <  r,  /=l,2,  r  =  0,  (5) 

where  a,  is  the  thermal  diffusivity,  is  the  thermal  conductiv¬ 
ity,  t  is  the  lime,  and  is  the  convection  heat  transfer  coefficient. 
The  above  boundary  value  problem  of  heat  conduction  is  solved 
by  applying  the  orthogonal-expansion  technique  over  multilayer 
regions  (Ozisik,  1968;  Tittle,  1965;  Tittle  and  Robenson,  1965). 
The  temperature  function  T(r,t)  can  be  represented  in  an  infi¬ 
nite  series  of  a  separated  space  (r)  and  time  (r)  variables  in 
the  following  form: 

r(r,  f)=  j=l,2,...,N  (6) 

n=l 

where  A„  are  arbitrary  constants  and  Xj„(r)  are  the  eigenfunc¬ 
tions  associated  with  the  above  eigenvalue  problem.  The  solu¬ 
tion  of  this  eigenvalue  problem  will  have  the  form 

Xj„(r)  =  BjnUoiujjnr)  +  Cj„Yo{uJj„r)] 

<  r  <  n  /  =  1,  2,  .  . . ,  (7) 


where  and  Y„  are  Bessel  function  of  the  first  and  second  kind, 
receptively.  Since  there  is  no  energy  storage  in  the  infinitesimal 
thickness  of  the  interface,  the  time  behavior  of  the  temperature 
at  the  interface  should  be  the  same  on  both  sides.  In  mathemati¬ 
cal  form,  this  condition  is  written  as 


=  a,ujj„  7  =  N  -  1,  .  .  .  ,  1.  (8) 


For  the  case  of  uniform  initial  temperature  (To)  all  over  the 
regions,  the  coefficient  is  written  as 


To*Xj,{r)^  r-dr 


To^Xjnir)-  r-dr 


(9) 


By  applying  the  boundary  conditions,  the  coefficients 
and  the  eigenvalues  cjy„  can  be  evaluated. 

Now,  consider  the  case  of  the  three  concentric  cylinders.  The 
temperature  functions  for  the  layers  are 


Ti(r,  r)  =  S  A„-X^n(r)-e  “■“■T-' 

1 

r2  <  r  <  r]  (10) 

TAr,  /)  =  I 

n=l 

r^  <  r  <  r2  (11) 

7-,(r,  f)  =  S  A„-X,Xr)-e 

0<r<r3.  (12) 

The  eigenfunctions  are 

X]„{r)  =  +  Ci,iFo(  u.'!„r)] 

^2  <  r  <  ri  (13) 

X2„ir)  =  B2n[Mi^'2.r)  +  C:M^'2.r)] 

ry  <  r  <  r..  (14) 

Journal  of  Applied  Mechanics 


JUNE  1998,  Vol.  65  /  347 


I 

1 

_ 

1 

t 

1 

1 

t 

_i _ 

I 

1 

Fig.  1  W-phase  concentric  cylinder  model 


materials  are  assumed  to  be  transversely  isotropic  both  in  stiff¬ 
ness  and  thermal  expansion,  where  (r,  is  the  transverse  plane. 

The  equilibrium  equations  in  the  r,  6,  and  z  directions  are 
given  by  (Boley  and  Weiner,  1960;  Manson,  1966;  Hasselman 
and  Heller,  1980;  Timoshenko  and  Goodier,  1971), 

dcJrr  1  r9  rz  ^ rr  ~  ^99  _  q 

dr  r  06  dz  ^ 


da«r 

1  da^e 

daoz 

+  ^  =  0 

dr 

r  de 

dz 

r 

da-r  1  BcTzo 

BcJtz 

+  ^  =  0. 

(23) 

dr 

r  06 

dz 

r 

The  stress-strain  relations  are 

(Trr  =  CwCrr  +  C[2^B9  +  GnC..  “  P\T 

(Tee  =  Cn^rr  +  CuCeo  +  -  /3\T 


0<r<r3  (15) 

The  eigenvalues  are  related  by 

C(\Lv]n  =  C(2i^2n  =  (16) 

By  applying  the  boundary  conditions  (Eqs.  (2),  (3),  (4)), 
one  can  get  the  following: 

—  h[jQ(ijO\nr\)  +  Ci„yo(^ln^*l  )]  (17) 

‘/o(t^’3«rO  =  B2n[JQ((^2nf’3)  +  C2„Tn(  )  ]  (18) 

^2rt[7o(^’2n^2)  G2nTo(^’2M^  2)] 

=  B}„[7o( Gi,|To(u.'ht^2)]  (19) 

k^UJ^J^iuJMjf'3)  =  /^2U^’2«^2/;[7l(u’2«rO  +  C:„Ti  (  ^'2«rO  ]  (70) 

k2U}2rtB2n['^l(^2nr2)  +  ( t^2/i''2 ) ] 

=  kiLu'[„B\n[J\{^\nf‘2)  +  CuKi  (  )  ]  .  (21) 

Equations  (17) -(20)  along  with  Eq.  ( 16)  could  be  solved  for 
the  following  unknowns,  C|„,  C^,,,  Bi,,^  u/S,,,  u,'?,,. 

The  solution  for  these  equations  is  obtained  by  eliminating 
all  of  the  unknowns  except  the  final  equation  is  solved  for 
the  root  cji,,  by  using  the  bisection  method.  Finally,  all  the  other 
unknowns  were  evaluated  in  terms  of  the  known  eigenvalue 
Lu’i,,.  The  coefficients  is  evaluated  by  using  Eq.  (9)  (;  =  3 
for  this  case).  After  integration  and  simplification.  A„  is  given 
bv 


(7-2  =  C^Crr  +  C\2^ee  +  6733€2: 

(T  rz  2G44C/-2 

(Jb:  =  2Cu^Bz 


(TrB  —  (6’n  “  C\2)^r0  (24) 

and 

I3\  =  (Cm  +  C!2)ar  +  C^a^ 

P2,  —  2C\2CiT  +  C330fi.  (25) 

where  Cij  are  elastic  constants,  and  aj  and  are  the  thermal 
expansion  coefficients  along  the  transverse  (r)  and  the  longitu¬ 
dinal  (z)  directions,  respectively.  The  strain-displacement  rela¬ 
tions  are  given  by 


_  1  /  d^  1  dll: 
"  2  r  dO 


1  \r 


(26) 

00  dr 


where  cr,^,  6,y,  and  w,  are  components  of  stress,  strain,  and  dis¬ 
placement,  respectively.  Due  to  the  axisymmetry,  the  displace¬ 
ment  field  can  be  expressed  as 


where  Z,(X„)  =  B„[J,(X,)  4-  CJ,{X„)],  Finally,  by  using  Eqs. 

{ 10)- ( 12)  the  temperature  field  is  well  defined. 

Stress  Field  in  N  Concentric  Cylinders  Subjected  to  Ther¬ 
mal  Loading.  Consider  the  N  concentric  long  circular  cylin¬ 
ders  shown  in  Fig.  1 .  The  outermost  cylinder  radius  is  r,  and  the 
innermost  cylinder  radius  is  i\  .  The  entire  model  is  subjected  to 
axisymmetric  temperature  change  ATi  r.  r ).  In  general,  all  the 


u^=i{{r).  0,  ii:^\v(z).  (27) 

A  substitution  of  Eq.  (27)  into  Eq.  (26)  yields  to 


It 


— 


r 


Dw 

Ik.' 


—  6b-  —  6,m  —  0. 


(28) 


Substituting  Eq.  (28)  into  Eq.  (24)  gives 
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^  du  u  dw 

(Jrr  “  ^11  ^12  “  ^13  “T  PM 

dr  r  oz 

c'  4-  r^  ^  4-  c  fi  T 

(Tm  =  (-12  TT  ^11  "  ^’3  ~  PU 

dr  r  dz 

9u  u  dw 

dr  r  dz 

Crz  =  (Tez  =  ^rB  =  0. 


(29) 


Then  the  governing  equations  are  obtained  in  terms  of  displace¬ 
ments  from  Eqs.  (23)  and  (29) 

C„^-C„-^^  +  ^^  =  /3,^,^  =  0.  (30) 


Now,  when  the  materials  are  isotropic,  both  in  stiffness  and  in 
thermal  expansion,  we  have 

C\\  =  C33  =  2/ii  +  \i,  CI2  =  CI3  =  X] 

Cli  =  C33  =  2^2  Xt,  C12  “  ^13  =  Xt 

(2ii  =  C33  =  2/^3  +  X3,  C12  =  C'n  =  X3 

P\  =  p\  =  {2r),  -h3X,)«i 

P\  =  I5\^  (27?2  +  3X2)«2 

(2t]j  +  3X3)a3-  (26) 

After  substituting  Eqs.  (32),  (34),  and  (36)  into  the  bound¬ 
ary  conditions.  Eq.  (31)  yields 


dr'’ 


r  dr  dr  dz 


Now,  for  the  case  of  three  concentric  long  circular  cylinders, 
the  boundary  conditions  are  given  by 

(y\r{r\)  =  0 

U\  -  M2,  vv^i  =  ^2.  CF\r  ~  crlr  at  r  =  r2 
M2  =  M3,  vvs  =  W3,  cr^r  =  at  r  =  r^, 

J  crl-jdr  -H  J  al-rdr  -!-  J  airdr  -  0.  (31) 

The  general  solution  to  Eq.  (30)  is  given  by 
M(r,  0  =  Ar  +  ^  ^  j  g(r,  t)dr  -  ^  j  r-g{r,  t)dr 

w(z)  =  Ez  +  F  (32) 


+  Ci(/-,,  t) 


Di(n.O 

■> 

r\ 


B, 


Air^  + - h  r^Ci(r2,  t)  + 

ri  '  ^2 


+  C!3£-y0|r,(r,,?)  =  O 
Di(ri,  t) 


=  +  ^  +  r.C4r2,  t)  + 

r2  '  r2 


.  B2  .  ^  ,  ...  D2{r,.t) 

A->r^  +  —  0  + - 

“  ‘  ''s  ''s 


where 


=  A3r3  +  r3C3(r3,  t)  + 


D4r„  t) 


f)  =  — -7- 
Cii  dr 


(33) 


F  =  0.  and  substitution  of  Eq.  (32)  into  Eq.  (29),  yields 


3  +  C:t(ry,  t)  — 


DAr,,  t)' 
n 


<7  rr  —  t-  I 


H-  Ch 


+  Cl 


/(j  +  C}{ry,  /) 


DAr>.  r) 


+  C„(r.  l)  + 


D„(r.  t) 


+  C]yE  -  l3\Tr{r„  t) 
A 


—  C'i2 


An  ~  — T  +  C„{t\  t)  — 


+  CliE  -  P'lTJr.  t) 
D„{r.  tf 


=  Ct> 

+  C-n 


B-,  ^  ,  D,(r,./)1 

3 - f  +  C.f/-,,  1) - ; - 

r?  ''5  J 

5.  „  ,  Z):(r<.  t)l 

Ai - f  +  £2(^3,?)  , 

r-y  r-,  J 


.  B„  ,  D„{r.t) 

+  Cn 

A„  H — 7  +  C„ii ,  r)  + 

Bz  ^  .  D.fr,,  ?)■ 

r-  fr- 

C]] 

A->  —  — 7  +  C'>(r-»,  /) 

_  ’  '*2  f'2  , 

+  C\iE  —  P'lTzii'j,  t) 


+  CUE  -  P‘\T„(r.  t) 

a",,  =  ICAAn  +  C„(r.  t)]  +  CUE  -  0'^T„(r.  t).  (34) 

where  n  is  referred  to  as  the  cylinder  under  consideration,  n  = 
1  is  the  outermost  cylinder,  n  =  2  is  the  middle  cylinder,  and 
n  =  3  is  the  innermost  cylinder,  and  where 


+  CU 


B2  .  ^  0:(':- 

A.  —  — ~  C'f{r->,  t)  ; 

rt  ■  ■  f; 


£13 £  “  P~\Fz{r2^  t) 

B, 


C„(r,/)4J 

D„{r.  f)  =  -^  J  r-gn(r. 


g„(r.  t)dr  =  T„(r.  t) 

1 1 

t)dr  =  L  A,„B„,„e-- 

-C  1 1  I 


Ddi'z.  t) 


fi,  ,  D,(r;.f) 

A| - r  +  Ci(r3.  t) - ^ - 


X  C,„„Y2(  /)]].  (35) 


+  Cl3£  -  P\Tdr2.  t) 
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Fig.  2  Dimension  of  TBC  layer  in  relation  to  bond  coat  and  superalloy 
substrate 


case  where  the  specimen  was  plasma-sprayed  in  air  with  a  thin 
zirconia-yttria  (Zr02-8  wt%  Y2O3)  layer  on  a  nickel-chro¬ 
mium-aluminum-zirconium  bond  coat,  as  in  Chang  et  al. 
(1987).  The  specimen  was  cooled  down  from  a  processing 
temperature  of  1000°C  to  a  temperature  of  0°C.  A  convection 
heat  transfer  coefficient  (h)  of  a  value  50  W/mVK  was  used 
in  the  analysis.  The  finite  element  model  was  used  for  the  case 
of  three  concentric  cylinders  using  ALGOR  for  a  fully  elastic 
analysis.  A  quarter  model  was  used  since  the  specimen  is  axi- 
symmetric  and  also  symmetry  is  maintained  about  the  middle 
section.  A  four-node  isoparametric  element  type  was  used;  the 
total  number  of  elements  was  500.  Along  the  longitudinal  axis, 
the  nodes  can  move  in  the  axial  direction  only,  and  along  the 
midsection  symmetry  line,  the  nodes  can  move  in  the  radial 
direction  only. 
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t) 
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Equation  (37)  has  six  linear  equations  with  six  unknowns 
{Al,  A2,  A^,  Bi,  B2,  E),  and  after  solving  these  equations,  the 
stress  field  is  fully  defined  at  every  r  and  r. 

Finite  Element  Stress  Analysis 

The  uncoated  nickel-based  superalloy  specimen  had  a  radius 
of  9.26  mm.  the  bond  coat  layer  thickness  is  0.184  mm,  the 
TBC  thickness  is  0.556  mm.  and  the  length  of  the  specimen  is 
100  mm.  These  dimensions  are  shown  in  Fig.  2.  The  properties 
of  these  three  layers  are  given  in  Table  1 .  We  analyzed  a  model 


Results 

The  orthogonal-expansion  technique  over  a  multilayer  region 
was  used  in  order  to  determine  the  temperature  distribution  in 
the  concentric  cylinders.  The  first  three  roots  (tjjn)  were  found 
to  be  dominant  in  prescribing  the  temperature  field.  All  other 
roots  are  weak  and  neglected.  The  axisymmetric  condition  was 
applied  on  the  equilibrium  equations  along  with  stress-strain 
relations  and  strain-displacement  relations,  so  the  equilibrium 
equations  reduce  to  two  ordinary  differential  equations.  The 
transient  stress  field  was  obtained  by  applying  the  boundary 
conditions  to  the  solution  obtained  from  those  ordinary  differen¬ 
tial  equations.  Transient  residual  stresses  have  been  evaluated 
for  the  case  of  three  concentric  cylinders  by  the  finite  element 
approach  as  well  as  by  using  the  elasticity  approach.  Figure  3 
shows  the  variation  of  stress-field  components  versus  time  at  r 
=  2  mm  in  the  substrate  region,  it  is  clear  that  the  stresses  are 
higher  for  the  case  of  a  partly  cooled  compared  to  a  fully  cooled 
case.  Because  of  the  high-temperature  gradient  in  the  super¬ 
alloy,  when  the  time  reaches  16.5  seconds,  the  stresses  are 
maximum.  Figure  4  shows  the  variation  of  stress-field  compo¬ 
nents  versus  time  at  r  =  9.3  mm  in  the  bond  coat  region.  It  is 
clear  that  the  stresses  are  higher  for  the  case  of  a  fully  cooled 
compared  to  that  of  a  partly  cooled  case.  Figure  5  shows  the 
variation  of  stress-field  components  versus  time  at  r  =  9.6  mm 
in  the  TBC  layer.  It  is  clear  that  the  stresses  are  higher  for  the 
case  of  a  fully  cooled  compared  to  a  partly  cooled  case  because 
the  temperature  field  has  a  high  gradient  at  the  earlier  stages 
compared  to  latest  stages.  The  temperature  distribution  as  a 
function  of  time  and  radial  distance  is  shown  in  Fig.  6.  It  is 
clear  that  the  temperature  field  has  a  high  gradient  at  the  earlier 
stages  compared  to  latest  stages. 

Figure  7  shows  the  variation  of  the  radial  stress  versus  the 
radial  distance  at  different  times.  The  peak  value  occurred  close 
to  the  center  at  r  =  16.5  seconds.  Also,  Fig.  7  shows  that  the 
radial  stresses  vanish  at  the  outer  surface.  Radial  stresses  are 
not  critical  for  this  specimen  configuration  because  of  their  low 
values.  Note  that  the  radial  stress  range  is  0-25  MPa.  Figure 


Table  1  Material  properties  at  22,  566  and  1149°C  (Wright,  1995) 


Material 

Young’s 

modulus 

(GPa) 

Poisson’s 

ratio 

Coefficient 
of  thermal 
expansion 

rc-') 

Coefficient 
of  thermal 
conductivity 
(W/m/k) 

Coefficient 
of  thermal 
diffusivity 
(m‘/s) 

Substrate 

175.8 

0.25 

13.91  X  !()-'' 

7 

2  X  10-" 

150.4 

0.2566 

15.36  X  lO-" 

13 

3.4  X  10-" 

94.1 

0.3224 

19.52  X  10"'' 

22 

2.64  X  10-" 

Bond  coat 

137.9 

0.27 

15.16  X  IQ-" 

25 

8  X  10-" 

121.4 

0.27 

15.37  X  10-'' 

50 

13.56  X  10-" 

93.8 

0.27 

17.48  X  10-" 

65 

15.52  X  10“" 

TBC 

27.6 

0.25 

10.01  X  10“" 

1 

5  X  10-’ 

6.9 

0.25 

11.01  X  10-" 

0.96 

3.51  X  10-' 

1.84 

0.25 

12.41  X  10" 

1 

3.43  X  10-’ 
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Fig.  3  Variation  of  stress-field  components  versus  time  in  the  superaltoy 

substrate  Fig.  5  Variation  of  stress-field  components  versus  time  in  the  TBC 


-  Radial  Stress 

- Tangential  Stress 

- Axial  Stress 


Time,  Sec. 


Fig.  4  Variation  of  stress-field  components  versus  time  in  the  bond  coat 


8  shows  the  variation  of  the  tangential  stress  versus  the  radial 
distance  at  different  times.  From  this  figure  one  can  see  that 
the  maximum  stress  occurred  in  the  bond  coat.  Figure  9  shows 
the  variation  of  the  axial  stress  versus  the  radial  distance  at 
different  times,  and  is  the  same  as  in  Fig.  8,  the  maximum  axial 
stress  occurred  in  the  bond  coat.  These  figures  show  excellent 
agreement  between  the  finite  element  solution  and  the  elasticity 
solution  for  the  steady- state  case.  Bond  coat /TBC  interface  is 
clearly  under  the  biaxial  stress  field.  The  highest  tensile  stress 
is  242  MPa  and  the  highest  compressive  stress  is  132  MPa. 
Figure  10  shows  the  variation  of  the  axial  stress  versus  the 
radial  distance  at  different  times.  In  this  figure,  temperature- 
dependent  properties  were  used  and  the  results  were  obtained 
by  using  the  finite  element  code  ALGOR.  By  comparing  this 
figure  with  Fig.  9,  it  is  clear  that  using  temperature-dependent 
properties,  will  increase  the  level  of  stresses  by  a  factor  of  1.7 
and  the  trend  will  remain  same.  This  effect  of  using  tempera- 
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Fig.  6  Temperature  versus  radial  distance  at  different  time  for  the  TBC/ 
bond  coat/ supera Hoy  system 


ture-dependent  properties  on  the  stresses  variation  for  all  of  the 
other  cases  were  the  same. 

Elasticity  solution  is  not  valid  close  to  the  free  edge  since 
the  boundary  conditions  are  not  satisfied  precisely  at  the  free 
edge.  The  boundary  conditions  are  subscribed  in  the  average 
form  at  any  cross  section,  which  may  not  be  precisely  valid  at 
the  free  edge. 

Compressive  stresses  in  the  TBC  play  a  major  role  in  TBC 
failure.  For  example,  a  flaw  size  of  12  times  the  TBC  thickness 
will  buckle  and  spallation  of  the  TBC  may  occur.  The  radial 
stress  across  the  interface  generates  a  driving  force  which  may 
initiate  an  interface  crack.  Subsequently,  the  interface  crack 
may  propagate  to  failure. 
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Fig.  7  Radial  stress  versus  radial  distance  at  different  times  for  the 
TBC/bond  coat/ superalloy  system 
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Fig.  8  Tangential  stress  versus  radial  distance  at  different  times  for  the 
TBC/bond  coat/superailoy  system 


Summary  and  Conclusions 

The  orthogonal-expansion  technique  was  utilized  over  a 
multilayer  region  to  determine  the  true  temperature  field  and 
subsequently  to  obtain  the  transient  residual  stress  field.  Analy¬ 
sis  of  the  concentric  cylinder  model  representing  the  TBC/bond 
coat/superailoy  system  showed  that  the  elasticity  solution  and 
numerical  results  from  finite  element  stress  analyses  (steady- 
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Fig.  9  Axial  stress  versus  radial  distance  at  different  time  for  the  TBC/ 
bond  coat/superailoy  system 
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Fig.  1 0  Axial  stress  versus  radial  distance  at  different  times  for  the  TBC  / 
bond  coat/superailoy  system  with  temperature-dependent  properties 


state  case  only)  correlate  well  in  predicting  thermal  residual 
stresses  in  the  various  layers  and  at  the  interfaces.  Because  of 
the  tree-edge  condition  and  stress  singularity  due  to  a  material 
mismatch,  the  elasticity  solution  is  not  valid  at  the  free  edge  of 
the  concentric  cylinders.  The  highest  stresses  were  found  to 
occur  at  the  bond  coat/TBC  interface.  At  this  interface,  results 
indicate  a  biaxial  stress  field  due  to  high  axial  and  tangential 
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stress;  the  radial  stress  component  was  relatively  small.  For 
transient  thermal  residual  analysis,  the  closed  form-elasticity 
solution  can  be  adequate  except  at  the  free  edge.  The  residual 
stresses  developed  in  the  substrate  are  higher  for  the  case  of  the 
partly  cooled  specimen  compared  to  the  fully  cooled  specimen 
because  the  temperature  field  has  a  high  gradient  at  the  earlier 
stages  compared  to  later  stages. 
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